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ABSTRACT 



We review the theory of polarized structure functions measured in deep- 
inelastic lepton-nucleon scattering, focusing on the most recent developments. We 
concentrate on the structure function g\, emphasizing the phcnomenological prob- 
lems related to the extraction of the proton and neutron gi{x,Q 2 ) from the data, 
and to the determination of its moments, especially the first moment. In particular, 
we discuss the theoretical and experimental uncertainties due to the Q 2 dependence 
of the data, small and large x extrapolations, QCD loop corrections, higher twist 
corrections, and, in connection to neutron experiments, nuclear effects. We also 
discuss the current status of the sum rules satisfied by the first moment of g\ and 
their theoretical interpretation. 



1. Introduction 

Polarized deep-inelastic scattering provides a handle on the matrix elements 
of spin-dependent operators. Because of the intricacies of spin physics, this poses 
exacting challenges to experimentalists and theorists alike. On the theoretical side, 
it turns out to be equally hard to extract meaningful physical information from the 
data, and to understand the significance of that information. 

This is perhaps best appreciated by recalling the present status of the mea- 
surement of the first moment of the polarized structure function g\\ 



Among polarized observables, this is the simplest to understand theoretically be- 
cause there exists only one operator with the appropriate spin at leading twist in 
the operator-product expansion, the flavor singlet fermionic axial current 



A measurement of Ti thus determines (in a way to be discussed extensively below) 
the (forward) matrix element of j§ in the target: 




(1.1) 




(1.2) 



{p,s\jg\p, s) = Ms'*AE, 



(1.3) 
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where p 11 , M, and s M are, respectively, the target four- momentum, mass, and spin 
(normalized as s^s^ = — 1). 

The experimental value of AE for a proton target, as given by the various ex- 
perimental collaborations which have measured it, is displayed in Fig. la. All the 
determinations are compatible within errors; however, the spread of the results in- 
dicates how hard these measurements are, and how carefully they should be taken. 
Actually, the experimental result which has been around longest 2 has been reanal- 
ysed by various theoretical groups: the results, shown in Fig. lb, again display a 
remarkable degree of variation (also in the error estimate). Finally, it is clear (albeit 
with hindsight) that whatever the measured value of AE, its interpretation will not 
be obvious: the current j§ is not conserved because of the axial anomaly, hence its 
matrix elements do not directly measure a well-defined conserved quantum number 
of a physical state. 

The increase in experimental precision calls now for a careful reanalysis of the 
arguments which go into the extraction of structure functions from the data: indeed, 
the theoretical uncertainty starts to be comparable to the experimental error. On 
the other hand, after the considerable amount of theoretical work which followed the 
1988 EMC experiment, 2 the subtleties involved in the interpretation of the results 
displayed in Fig. 1 are now largely understood. 

In this paper we will concentrate on the phenomenological aspects of structure 
function measurements, with particular regard to recent developments, and we will 
briefly summarize the present understanding of their theoretical import. In Sect. 2 
we will provide some background on the polarized structure functions gi(x, Q 2 ) and 
g2(x, Q 2 ). In Sect. 3 we will then proceed to a discussion of the phenomenology of 
the proton structure function g\: first we will discuss the way g\ is extracted from 
the data, and review the problems involved in Q 2 dependence and extrapolation in x; 
then we will examine the way AE [Eq. (1.3)] is evaluated once gi is known, discussing 
the dependence on weak decay constants and the role of QCD loop corrections. In 
Sect. 4 we will review the current status of the theoretical interpretation of the data 
on AE, and we will discuss parametrization of polarized parton distributions. In 
Sect. 5 we will review problems which are of special relevance for the measurement 
of the neutron structure function g™, in particular nuclear effects and higher twist 
effects, and in Sect. 6 we will discuss the significance of the proton and neutron 
measurements taken together.* 

2. The structure functions and their moments 

Polarized structure functions are the form factors which parametrize the cross 
section spin asymmetry for deep-inelastic scattering of polarized leptons off a po- 

* For a more detailed introduction to the theory of g\ in QCD and to theoretical controversies 
spurred by its measurement the reader is referred to earlier reviews, 7,10 while a recent assessment 
of the theoretical status of the g\ measurement is in Ref. 11; for a detailed review of the theory 
of gi see Ref. 12, and for a discussion of other polarized structure functions (relevant for different 
processes, such as Drell-Yan scattering) Ref. 13. 
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larized hadronic target 
hadronic tensor 15 ' 14 



d ^dQ^dv — they are given by the antisymmetric part of the 



iW%" = Jd 4 xe iq - x (p,s\J^(x)J^(0)\p, s) 



%Me vvpa q p 



—9i(x,Q 2 ) + 



p-q 



( P - q y 



(2.1) 



(with standard parton model kinematics), where M is the target mass. Polarized 
neutrino scattering is only of academic interest at present, hence we will neglect 
weak interaction effects, and assume J M in Eq. (2.1) to be electric currents. 

The light-cone expansion of the current product in Eq. (2.1) implies that the 
moments of the structure functions g± and gi i n the Bjorken limit are given in the 
nonsinglet case by 15 ' 14 



/ 

Jo 

f 

Jo 



dx x 



n-l 



gi (x,Q 2 ) = -C?(Q 2 )a n (Q 2 ) 



(2.2) 



dxx^g^Q 2 ) = IL— [C2(Q 2 )d n (Q 2 )-CnQ 2 )a n (Q 2 )} . (2.3) 



C n (Q 2 ) are perturbatively calculable coefficient functions, and a n and d n are given 
by matrix elements of the leading twist polarized operators: 

Ma n s {a p^ ...p" n - l} = -{p^li^^^D^ ...D^- l} X^\ps} (2.4) 
Md n s {a p^ ...p" n - l} = (p,s|i n - 1 ^757 {[fJ J D m] ...D^- l} \iiP\ps), (2.5) 

where {} denotes complete symmetrization. Notice that, even though the operators 
Eq. (2.2) are twist- 2 and the operators Eq. (2.3) are twist-3, their respective con- 
tributions to the light-cone expansion of W^ v are of the same order in Q 2 . In the 
singlet case, these operators will further mix with gluonic ones;* 16 however, the first 
moment of gi, which we will be mostly concerned with, is still given by Eq. (2.4). 

The reason why we will concentrate our discussion on gi is that for longitudi- 
nally polarized protons the only nonvanishing components of the hadronic tensor 
are 14 (in the frame where, say, only the third spatial component of q is nonzero and 
its energy vanishes) 



W 



12 

A 



-w 



21 
A 



± 



9i ~ 



(2Mx\ 



92 



(2.6) 



* Note also that Eq. (2.5) is true only for strictly massless quarks; in general the operators con- 
tributing to 32 will mix with twist-3 mass dependent operators. The contribution of these operators 
does not vanish asymptotically, even though their matrix elements should be of order m/M, with 
m and M the quark and nucleon masses. 
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hence only g\ is relevant asymptotically. For transverse polarization (along, say, 
the first spatial coordinate) the nonvanishing components are instead 



Wf = -Wf = 



t20 



2Mx 



[91 + 92] , 



(2.7) 



thus, even though the two structure functions contribute equally to it, the whole 
cross section vanishes asymptotically. Only an experimental upper bound on gi is 
thus presently available; 17 an alternative determination of it is forthcoming. 6 

In general, Y\ will receive both singlet and nonsinglet contributions, and it is 
therefore convenient to define axial charges a* for the i-th flavor: 



Mcus^ = (p, slipa^ipilp, s). 
The first moment of g\ [Eq. (1.1)] for a nucleon target is then given by 



(2.8) 



-nP, n 
1 1 



(Q 2 ) 
1 

12 



1 

2 



N f 



i=l 

2' 



C NS {Q 2 ) 



± (a u - a d ) + - (a u + a d 



2a, 



+ C S (Q 2 )- (a u + a d + a s 



(2.9) 

where in the last step we have assumed that only the three lightest flavors are 
activated, the plus (minus) sign refers to the proton (neutron), and the explicit 
form of the singlet and nonsinglet coefficient functions Cs, Cns will be given in 
Sect. 3.2. 



3. Phenomenology of g\ 

The primary quantity from which experimental information on polarized struc- 

T T T J- 

ture functions is obtained is the polarization asymmetry A(x,Q 2 ) = ^ TT ~^ n . In 
order to relate this to W^ v [Eq. (2.1)], we must extract from A the asymmetry for 
scattering of transverse {i.e. with longitudinal helicity) virtual photons, defined as 

M x ,Q>) = ai < 2 - a: >< 2 (3.1) 

01/2 + 03/2 

where the subscripts denote the total angular momentum of the photon-nucleon 
pair along the incoming electron's direction; the denominator of (3.1) is the total 
transverse photoabsorption cross section ot = 0"i/2+°3/2- A is related to A\ by A = 
D(Ai+r]A2), in terms of the factor D which provides the longitudinal depolarization 
of the virtual photon with respect to its parent lepton, the interference A2 with the 
longitudinal photon polarization amplitude, and a kinematic factor 77. 
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The asymmetry A\ measures the ratio between the combination (2.1) of po- 
larized structure functions, and the unpolarized structure function F\ (the other 
independent unpolarized structure function gives the total longitudinal cross sec- 
tion). In the Bjorken limit the contribution of (72 is negligible [Eq. (2.6)] only g\ 
contributes to the numerator in Eq. (3.1), and g\ is given by 

g 1 (x,Q 2 ) = A 1 (x,Q 2 )F 1 (x,Q 2 ) (3.2) 

(3 - 3) 

where the last step holds at leading twist. The unpolarized structure function F2 
is directly determined from experiment, and has the simple parton model interpre- 
tation (which in the DIS scheme remains true to all orders in perturbative QCD) 

-^- = Y,e 2 [q i (x) + q i (x)}, (3.4) 
i=l 

while R = — , where ol is the longitudinal virtual photoabsorption cross section. 
Because in the Bjorken limit R vanishes, Eq. (3.4) reduces to g\ ~ A\Fil^2,x)\ in 
particular, in the naive parton model it thus follows that 

<7r rt °» = E e ' [Agi(x) + Aft(x)] , (3.5) 

i=l 

where Aqi are the polarized parton distributions for quarks of flavor %. 

A precise determination of g\ thus hinges on the following assumptions. First, 
A\ must be determined from the measured asymmetry A: in older experiments 1 ' 2 
this was done by simply neglecting A2 and identifying A with DA\. Then, the 
contribution of gi to A\ is neglected, exploiting Eq. (2.6), so that Eq. (3.2) holds. 
Finally, g\ is determined from the measured F2 using Eq. (3.3), the main uncer- 
tainty being in the experimental knowledge of R. The first assumption is under 
control: A2 satisfies a positivity limit j^l < v^R, and is now 17 measured to be 
compatible with zero within uncertainties, so that (taking also advantage of the 
fact that the coefficient rj turns out to be small in the kinematic range covered by 
present experiments) its possible effects are just included in the systematic error. 
The neglect of (72 should not affect the determination of the first moment of gi, 
since, regardless of the size of g2, its first moment is expected to vanish (according 
to the Burkhardt-Cottingham sum rule); 12 ' 18 even if it did not (due to nonpertur- 
bative effects) its contribution is estimated 19 to modify the first moment of g\ by at 
most 1%. A precise measurement of (72 would settle the issue. The last step of the 
analysis is also under complete control: due to cancellations in the dependence of 
D and F\ on R, the asymmetry gi is in fact essentially independent of R. 4 All these 
effects combined lead thus to an uncertainty of a few per cent, which is included 
in the experimental error and is smaller than the uncertainty on the unpolarized 
structure function F 2 . 

In the remainder of this section, we will consider in particular the proton ex- 
periments, which afford the best experimental accuracy. 
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3.1. Scale dependence and x- extrapolation 

The values of g\ which can be determined using Eq. (3.3) span in practice a 
limited kinematic range: very small and large x regions are excluded, while low x 
data are taken at low Q 2 and conversely. Thus, in order to obtain a determination 
of gi(x,Q 2 ) for all x and one single value of Q 2 , the data must be evolved to a 
common Q 2 and then extrapolated over the whole x range. 

The scale dependence of the asymmetry A\ is given by the Altarelli-Parisi 
equations, and follows from the fact that g\ and depend on Q 2 in different ways, 
even at leading order. This is apparent if one considers the respective first moments: 
at one loop the anomalous dimension of Ti vanishes (the axial current is classically 
conserved; this conservation is spoilt by the axial anomaly, but this is a two-loop 
effect), while the first moment of F2 evolves since obviously the total number of q-q 
pairs is not conserved. The scale dependence is thus calculable, the only theoretical 
uncertainty being that related to the relatively poor knowledge of polarized parton 
distributions (see Sect. 4.2). The result of the calculation, 20 with two extreme 
scenarios for the polarized gluon distribution, is displayed and compared with the 
data of Ref. 1 in Fig. 2a. 

Clearly, the data are not precise enough to display the calculated scale depen- 
dence, even though they are consistent with it (the same applies to more recent 
data 3 ), which explains why the scale dependence of the asymmetry is not yet ob- 
served directly. Nevertheless, the average scale of the smallest and largest x bins is 
often rather different: for example in Ref. 4 the smallest x bin ((x) = 0.005) has 
(Q 2 ) = 1.3 GeV 2 , and the largest x bin ((x) = 0.48) has (Q 2 ) = 58.0 GeV 2 . It is 
therefore important to correct the data for this effect, i.e. to evolve the measured 
asymmetry to the common scale at which gi is to be determined. The size of the 
correction (shown in Fig. 2b for the data of Ref. 2) depends strongly on this scale: 
the correction on is smaller than 5% for the determination of Ref. 2, which is 
given at Q 2 = 10.7 GeV 2 , but is for instance of order 10% for the deuteron data 22 
at Q 2 = 4.7 GeV 2 . The overall effect is anyhow larger than the uncertainty in- 
volved in its determination (which is mostly due to the uncertainty in the form of 
the polarized gluon distribution) and should therefore be included in the analysis 
of the data, especially if g\ is to be determined at low (Q 2 < 5 GeV 2 ) scales. 

Once the data are evolved to a common scale, g\ is still only known over a 
limited x range: for example, the combined data of Refs.21 and 2 have 0.01 < x < 
0.7, and the more recent ones of Ref. 4 have 0.003 < x < 0.7. In order to be able 
to compute moments of g±, which are the primary quantities in QCD, gi thus has 
to be extrapolated at small and large x. 

The large x behavior should be controlled by QCD counting rules: 23 ' 24 because 
for large x the virtuality of the struck quark is large, the scattering process should 
be purely perturbative, and the dominant contribution to it should come from 
minimally connected tree graphs, i.e. those where three valence quarks exchange 
two perturbative gluons. This leads to the prediction 

G q ~(l- X ) p ; p = 2n-l + 2As z , (3.6) 
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where G q is any parton distribution, n is the number of spectator quark lines (i.e., 
n = 2 for a nucleon) and As z — (As z = 1) for parallel (antiparallel) quark and 
proton helicities. This implies that in the approximation of Eq. (3.5) one would 
expect 

A x ~ const.; g x ~ (x - 1) 3 . (3.7) 

x — >1 x — >1 

The fall-off of g\ at large x implies that the contribution of the large x ex- 
trapolation to Ti is expected to be small, i.e. around a few per cent of T\. The 
extrapolation has been performed by fitting a phenomenological function either to 
Ai (which amounts to fitting a straight line through the last few data points), 2 or 
to <7i. 5 The (preliminary) result of the latter extrapolation actually turns out to be 
surpisingly large (by a factor 3 larger than that of Ref. 4) ; this appears to be due to 
an actual difference in the data, and not to the method used in the extrapolation; 
however, the large x points have large uncertainties — indeed all experiments are 
compatible with each other — and the effect on Ti is anyway modest. The same 
applies to results obtained 24 by imposing the behavior Eq. (3.6). A different esti- 
mate has also been obtained by using a valence quark model, 9 which is consistent 
with the behavior (3.6) but predicts also its normalization; the results are similar 
(see however the discussion of the deuterium data in Sect. 5). 

The small- a; extrapolation is somewhat harder to control. Arguments based on 
the dominance of known Regge poles 25 lead to the expectation 

9i ~ x a ; 0<a< 0.5. (3.8) 

x— »0 

However, use of perturbative QCD together with an Ansatz for the nucleon wave 
function 24 suggests that the gluon distribution should behave as 

^ ~ *i (3.9) 

(jr X— >0 

thus, if the unpolarized gluon distribution is dominated by a soft pomeron, 26 then 
G(x) ~ - so that the lower bound for a in Eq. (3.8) is saturated; but if G(x) 
has a harder behavior (such as a supercritical or perturbative pomeron 26 ) then 
a < (a ~ —0.1 with a supercritical pomeron, and even lower with the Lipatov 
hard pomeron 26 ). Due to QCD evolution the gluon feeds into the quark at small x 
(the eigenvector of the evolution being a linear combination of Aq and Ag), and a 
singular behavior with the same value of a is then induced in g\. 27 Also, a model of 
the pomeron based on nonperturbative gluon exchange 28 leads to the still singular 
but softer behavior 

g 1 ~ -2 lux. (3.10) 

x^O 

Finally, it has been argued 29 that negative signature cuts could induce an even more 
singular behavior 

9i ~ -rV, (3-11) 

X1I1 X 
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even though there is no compelling theoretical evidence in favor of such contribu- 
tions. It is worth noticing that not all these behaviors are stable upon perturbative 
evolution: 27 in general if, at some scale Qo, gi is regular, or has a singularity of the 
form (3.8) with a < 0.1 or softer [such as (3.10)], then at larger Q and small x, due 
to QCD evolution, it will develop a double logarithmic singularity 30 

gi ~ o exp2 7v /logl/x log log Q 2 ; 7 2 = (16/33-2ra/) [5+4^1 - 3n//32]. (3.12) 

More singular behaviors such as (3.11) will instead be preserved by the evolution. 
The issue is further complicated by the mixing of quark and gluon contributions, 
which turn out to contribute to the small-x eigenvector with opposite signs. 27 

Keeping in mind that the theoretical picture is somewhat blurred, let us turn 
to the data. The earlier experiment, 1 which extends down to x = 0.015, fitted to the 
asymmetry a parametrization that behaves at small x as A± ~ x 1 which then, 

x— >0 

assuming F 2 ~ const., leads, by Eq. (3.3), to the behavior (3.8) with a = 0.12. 

x— >0 

A direct fit of the form (3.8) to the last seven data points 31 leads to a = 0.07lo' 32 ,. 
Finally, a recent determination 32 of F 2 displays a fitted behavior F 2 ~ x~ , 

a;— >0 

which, if used together with the above parametrization of A±, would lead to 9 a = 
0.02; however if g\ is determined 3 from the measured asymmetry 1 using this form 
of F 2l the result can be still fitted with the function (3.8). Quantitatively, these 
methods all lead to a contribution to the first moment from the x < 0.01 region 



r s r 



. 0.01 

small x 



/•u.ui 

= / dxgx{x) « 0.002, (3.13) 
Jo 



which is about 2% of Y\ (and usually just included in the systematic error on the 
value of Ti). 

The more recent experimental data, 4 however, extend to significantly smaller x 
(xmin = 0.003), and lead to the striking result displayed in Fig. 3: not only there is 
no evidence of a fall-off with x, but actually (even though with large errors) the last 
four data points display a rise at small x. These data have been extrapolated by 
fitting a constant value for gi(x) to the last two points; the resulting contribution 
from the unmeasured small x region to Ti is now larger than in Ref. 2 by a factor 
2, and adding to that the contribution from the (now measured) 0.003 < x < 0.01 
range leads to a value of rf mallx larger than Eq. (3.13) by a factor 6 or 7. Fitting 
the behavior (3.10) leads 29 to a similar result, while the very singular form of 
Eq. (3.11) gives a value which is yet larger by a factor 3, i.e., about 20 times larger 
than Eq. (3.13). This alone would increase 29 the value of Fi by almost 25%. 

Of course such conclusions should be taken with caution, especially since what 
is actually observed experimentally is an approximately constant behavior of the 
asymmetry A±, and the rise is then produced by extrapolating a parametrization to 
F 2 outside its declared 32 range of validity; furthermore, perturbative evolution may 
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substantially affect the small x behavior 27 and it is unclear that a small x tail fitted 
at (Q 2 ) ~ 1 GeV 2 can be assumed to be unchanged at (Q 2 ) = 10.7 GeV 2 where gi(x) 
is given. Unfortunately the very recent SLAC data, 5 which have better statistics, 
only extend down to x m i n = 0.029; they are consistent with the EMC/SMC data 
and, if fitted to a constant, 6 lead to a value of rf malla: of the same size as that of 
Eq. (3.13). 

3.2. Determining AE 

The singlet and nonsinglet components of Ti can be extracted using Eq. (2.9). 
This requires the determination of the matrix element of the triplet and octet cur- 
rent, and the computation of the coefficient functions C(Q 2 ). 

The triplet and octet current are conserved and therefore scale independent; 
they can thus be taken from any other process. The triplet matrix element is 
accurately known, because (through trivial isospin algebra) it is equal to the axial 
coupling measured in nucleon /3-decay: 

a u -a d = ^- = F + D = 1.2573 ± 0.0028, (3.14) 

where F and D are octet meson decay constants in the standard SU(3) parametriza- 
tion. The octet matrix element is then given, using SU(3) symmetry, by 

a u + a d -2a s = 3F - D = 0.579 ± 0.025. (3.15) 

The value in Eq. (3.15) is obtained from the most recent fit 29 to octet decays; notice 
that these values are significantly different (and more precise) than those used in 
older analyses. 2 

Now, AE is found rewriting Eq. (2.9) as 

C S (Q 2 )AE(Q 2 ) = 9TI - Cns ^ ( 3 F + D) . (3.16) 

The two terms on the r.h.s. of Eq. (3.16) are roughly of the same size (typically, 
Ti ~ 0.1), thus AE arises from a large cancellation between them. In the most 
recent and accurate data 4 the error on Y\ is of order 10%, hence it still dominates 
over the error on 3F + D : which is of order 1%. 

Knowledge of the perturbative coefficient functions in Eq. (3.16) has consider- 
ably improved recently. The one-loop results have long been known, in both the 
nonsinglet 15 and singlet 16 cases; the nonsinglet two loop 33 and three loop 34 co- 
efficients have been determined subsequently; the singlet two loop coefficient was 
recently calculated, 35 ' 36 and estimates for the next corrections (i.e. three loop in 
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the singlet case, and four loop in the nonsinglet) have also been proposed. 37 The 
full results (in the MS scheme) are 



C NS (Q 2 ) = [l - (^) 



- 3.5833 




-20.2153 




(~ 130) (*) 



+ ... 



(3.17) 



Cs(Q>)= i-(^) 



- 1.0959 



+ ... 



(3.18) 



where ~ denotes the estimated coefficients. The estimates are arrived at by requir- 
ing minimal scale sensitivity (plus, in the nonsinglet case, some guesswork), and 
should be taken with care. Two comments are in order here: first, higher loop 
corrections turn out to have relatively large coefficients and are not negligible if 
Q 2 < 5 GeV 2 ; also, all loop corrections go in the same direction, namely, for given 
]?! the value of AE obtained from Eq. (3.16) increases as more perturbative or- 
ders are included. The overall effect of these corrections is of order 20% around 
Q 2 ~ 5 GeV 2 ; this is large enough for the uncertainty on Aqcd to reflect on an un- 
certainty on AE of a few per cent. This sensitivity can be actually used to measure 
a s (see Sect. 6). 

The scale dependence of the quantity on the l.h.s. of Eq. (3.16) is still not 
entirely specified, because the singlet matrix element AE depends on Q 2 (due 
to the anomalous nonconservation of the singlet axial current). Its scale depen- 
dence starts at next to leading order; the first two nontrivial coefficients have been 
computed, 16 ' 36 while the next order is estimated 37 : 



4. The proton experimental results and their meaning 

The current experimental knowledge on Ti and AE for the proton is summa- 
rized in Table 1; the last column gives the proton matrix element of the strange 
axial current, and thus provides a measure of the Zweig rule violation in this chan- 
nel. The first two rows of the column are obtained from the same raw data; the 
more recent value 3 of AE differs essentially because of the use of updated values 
of F and D and the inclusion of higher loop corrections (the latter, however, have 
a negligible effect at this scale). The values of AE shown in the table are as given 
by the respective references, except the value for the SLAC experiment, which we 
determined using Eqs. (3.16) with the coefficient functions (3.17)-(3.18) and the 
decay constants (3.14)-(3.15). As more data become available, however, it is now 



AE(Q 2 ) = fl + | f— ) + 1-2130 



3 V 7T J 



(v) 2 + <~->(v) 3 



+ ... AE(oo). (3.19) 
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Ref. 


m 


r? 


AE 


a s 


2, 21 


10.7 


0.126 ±0.018 


0.12 ±0.17 


-0.19 ±0.06 


3 


10.7 


0.126 ±0.018 


0.14 ±0.17 


-0.15 ±0.06 


4, 21 


10.0 


0.142 ±0.014 


0.27 ±0.13 


-0.10 ±0.05 


5 


3.0 


0.129 ±0.011 


0.28 ±0.11 


-0.10 ±0.04 



Tabic 1: Summary of proton experimental results. All results hold at the given scale. 

more sensible to perform a global fit of AE from the determinations of Y\ in various 
experiments; 11 we will discuss this in Sect. 6. 

All data points turn out to agree within errors; the SMC result 4 is substantially 
larger because of the unexpectedly large values of g± in the small x region, which is 
not covered by other experiments. The recent, more precise (but still preliminary) 
SLAC data 5 also agree with previous data; the difference in value of Ti is mostly 
accounted for by Q 2 evolution. The error on AE is dominated by the uncertainty on 
Ti; this, in turn, comes in roughly equal proportions from statistics and systematics 
in the EMC/SMC experiments, whereas it is systematics-dominated in the E143 
experiment, 5 which has better statistics but a smaller kinematic coverage. If the 
large values of g± at small x of Ref. 4 are confirmed, then the small-x extrapolation 
of Ref. 5 will have to be corrected, since the effect on F\ is then comparable to the 
total uncertainty. 

4-1. Theoretical interpretation 

The large violation of the Zweig rule displayed by the data of Table 1 has been 
variously dubbed "spin crisis" or "spin puzzle" . Why this result should be puzzling 
at all (after all the Zweig rule is only a phenomenological expectation, known to fail 
in some channels) is perhaps best understood from the point of view of someone 
who wishes to construct a parametrization of the polarized quark distributions Aqi 
and gluon distribution AG. In leading order, AE is scale independent; also, in the 
parton model [Eq. (3.5)] AE is just the fraction of the nucleon helicity carried by 
quarks. Parton model results are modified by QCD evolution, but since the first 
moment does not evolve this identification is retained. This is in keeping with the 
observation that [Eq. (1.3)] AE is the coefficient of proportionality between the 
nucleon helicity and the quark axial charge, which coincides with the quark helicity 
for free, massless quarks. Thus, AE is the normalization of the polarized parton 
distribution which is input to the perturbative evolution, and one would guess that 
it should just be equal to the spin fraction carried by quarks in quark models of 
the nucleon, which is typically 7 ~ 0.6 ± 0.1; hence the puzzle. Of course, this is a 
puzzle for parton models and parametrizations, but not for QCD, even less so for 
effective models of the nucleon. 

The resolution of this difficulty in the parton model 38 is easy to state, and 
becomes apparent when the QCD evolution equations are solved in next to leading 
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order. The eigenvector of the QCD evolution equation for the first moment turns 
out to be 

AE = E / A ^ = ^l dx - 2^ A ^) (41) 

whose anomalous dimension is equal to that of the axial current Eq. (3.19) (the 
eigenvector remains Eq. (4.1) even at higher orders). Now, the crucial point is 
that this anomalous dimension starts at next to leading order, which means that at 
leading order the combination f^Ag is scale invariant; thus, the gluon contribution 
to Eq. (4.1) is not asymptotically suppressed by a power of a s , and the next to 
leading order analysis tells us that the one loop scale invariant combination which 
in the parton model is associated to AE is that given in Eq. (4.1), and not the 
naive parton model one of Eq. (3.5). The polarized gluon distribution whose first 
moment appears in Eq. (4.1) is uniquely defined and can be independently measured 
in different hard processes, 39 ' even though in polarized deep-inelastic scattering 
it is only the combination in Eq. (4.1) which is measurable. It appears that the 
coefficient of the gluon contribution to Eq. (4.1) can be changed by using different 
regularizations of infrared collinear singularities, 41 thus modifying the amount of 
mixing between quarks and gluons, but a detailed analysis 42 shows that this can 
only be done by including soft contributions in a hard coefficient function. This 
said, the gluon coefficient in Eq. (4.1) can still be modified by reabsorbing the 
gluon contribution in a redefinition of the polarized quark distribution, i.e. by a 
change of scheme; however, the definition Eq. (4.1) is the unique one where Yli A (ft 
is conserved to all orders, as the quark helicity must be. 

The reason for these results resides in the anomalous conservation law satisfied 
by the axial current 

d l Jg = N f -±tTe>»""F lu ,F pat (4.2) 
Ztv 

which is exact nonperturbatively. Due to the anomaly, matrix elements of the axial 
current do not measure directly the quark helicity, but rather a combination of 
the quark helicity and an anomalous non-conserved contribution; the latter reduces 
perturbatively to the expression Eq. (4.1), but can also receive a nonperturbative 
contribution O, which need not have a partonic interpretation, so that in general 43 

AE = ^2 fdx Aqi - N f Q. (4.3) 

The anomalous contribution is due to the helicity generated by anomalous particle 
creation, 43 in analogy 44 to electroweak baryon number generation.* In principle, 

* The fact that the gluon contribution can be absorbed in the quark one by a choice of scheme, 
but at the expense of losing the conservation of first moment of the polarized quark density, cor- 
responds to the fact that the anomalous dimension of the axial current can be set to zero by a 
finite renormalization (because it starts at two loops), but the mixing of the axial currents with the 
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the perturbative anomalous contribution Ag can be measured directly 10 ' 11 , while 
the nonperturbative one O can be revealed in different processes; 46 the former is 
scale dependent, as discussed above, while the latter is not. 

Regardless of whether it is Ag or a nonperturbative contribution O, a large 
cancellation between the quark contribution and the anomalous gluon contribution 
to AE must be invoked to explain its observed small value. This cancellation in turn 
can only be explained in terms of a nonperturbative mechanism. While effective 
models of the nucleon, such as the Skyrme, 47 ' 48 bag, 49 or Nambu-Jona-Lasinio 50 
models can easily be made to accommodate the observed values of AE, F and D 
(but at the expense of intoducing extra free parameters!, ), in these models it is 
impossible to separate quark and gluon contributions to observed quantities. On 
the other hand, a cancellation between Aq and O in Eq. (4.3) can be explicitly 
shown to take place due to instanton-like configurations in the QCD vacuum, 51 
however it is hard to go beyond model calculations. 

A rather different understanding of the small value of AE is arrived at in 
t-channel approaches, where the nucleon matrix element of the axial current is 
computed coupling the current to composite operators that correspond to physical 
bound states, which then couple irreducibly to the nucleon. The small value of Ti 
is then due to a strong nonperturbatively induced scale dependence, 52 or to the 
smallness of the relevant bound state propagator, which can be expressed using an 
exact Ward identity in terms of the topological susceptibility of the vacuum, 53 and 
estimated using QCD sum rules, 54 with results in excellent agreement with exper- 
iment. These approaches, which rely only on the nonperturbative QCD dynamics, 
imply that the smallness of Ti depends on the structure of the axial current, and 
is unrelated to the specific target which is being considered; 54 their relation to the 
parton approach, however, is not immediate. 

4-2. Polarized parton parametrizations 

The problem of constructing a parametrization of polarized parton distributions 
is actually not an academic one. Earlier proposals, based on the idea of relating 
polarized distributions to unpolarized ones through a dilution factor 55 lead neces- 
sarily to results in disagreement 1 with the data; however, once the anomalous gluon 
contribution Eq. (4.1) is taken into account, the data can be fitted phenomenologi- 
cally without giving up the idea that polarized and unpolarized quark distributions 

(gluonic) anomaly operator cannot be removed (being a one- loop effect), and the coefficient of the 
mixing (the strength of the anomaly) is fixed in a scheme independent way once the way in which 
Aq renormalizes is specified. Of course, gluons will also contribute to higher moments of gi, but 
this contribution is entirely scheme-dependent; the pertinent coefficient functions are given within 
specific choices of scheme in Ref. 45. 

t For instance, the pure Skyrme model in the Nf — > 00 limit predicts 47 AS = and cannot fit 
F and D accurately, but can be made to fit all the data by introducing corrections in j^j such as 

those due to coupling to vector mesons. 48 
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should have the same gross features, either within the same approach, 6 or by con- 
structing a parametrization from scratch. 57 Of course the data can also be described 
without anomalous contribution but assuming a large amount of Zweig rule viola- 
tion (which is absent in the unpolarized case). 

As more data become available, however, it is now possible to construct po- 
larized parton distributions which incorporate all the constraints from perturbative 
QCD [Eqs. (3.6), (3.9)] as well as the general constraint from Regge theory that 
the small x behavior be isosinglet. 45 Imposing the constraints at a starting scale of 
Q 2 = 4 GeV 2 the data can then be fitted rather accurately, consistently including 
the scale dependence at leading order, and imposing cuts to keep higher twist cor- 
rections under control. The fit can be performed assuming the polarized sea quark 
distributions to be entirely generated by perturbative evolution; AG turns out to 
be largely unconstrained, and can be chosen to have different large x behaviors. If 
the gluon is assumed to behave as a constant at small x, then the fit to the proton 
data 21 ' 2 ' 4 favors a singular quark of the form Eq. (3.8) (which dominates g\ at small 
x), with a ~ —0.55 ± 0.15. The neutron structure function can then be predicted, 
the data not being good enough to constrain it seriously. Eventually, it should be 
possible to test these parametrization directly: while AG can be measured typically 
by photon-gluon fusion, 10 Aq can in principle be measured for sea and valence fla- 
vors separately, by tagging final state hadrons; preliminary results 58 are consistent 
with the results of the fit to the proton data. 45 



5. Polarized scattering on neutrons 

A measurement of g™ provides the determination of an independent linear com- 
bination of polarized quark distributions. Thus, neutron experiments, besides giving 
independent theoretical information on AE, also offer the possibility of measuring 
directly the nonsinglet polarized structure function g\ — gf, whose first moment is 
just [from Eq. (2.9)] 

r?-r? = ^ s (Q 2 )(a.-a d ). (5.1) 

A comparison of the data with the known value of a u — ad [Eq. (3.14)] and the 
computed Cns [Eq. (3.17)] thus allows a direct test of isospin in this channel, 
as well as of the predicted scale dependence. However, neutron experiments are 
somewhat subtler, both because nuclear effects have now to be taken into account, 
and because they generally have a more restricted kinematic coverage than proton 
experiments. 

The neutron structure function can only be measured by scattering on nuclear 
targets. Specifically, deuterium 22 ' 5 and 3 He targets 59 have been used. In the former 
case, assuming additivity, one obtains a determination of g\ + g™, in the latter case, 
the two protons are mostly in an S'-state so that the spin of the nucleus is carried 
by the neutron. In both cases, however, there are complications due to nuclear 
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structure. Taking into account that the deuteron can be in D wave with probability 
ujd one has 22 

The 3 He wave function has in general 10 components; a simplified description in 
terms of a three-component wave function leads to 60 

17 = (1.15 ± 0.02)rf e + (0.057 ± 0.009)r?, (5.3) 

which however only holds in the Bjorken limit, and could be significantly corrected 
at finite Q 2 ; the error involved is still negligible compared to present day experi- 
mental accuracy. Further nuclear effects are due to Fermi motion (which for the 
deuteron is estimated 61 to lead to a sizable correction to gi, up to perhaps 10%), 
and to shadowing and antishadowing. 62 No systematic investigation of the latter 
effects is available yet. 

The problems related to Q 2 evolution and x extrapolation discussed in Sect. 3.1 
are even more serious in neutron experiments. An example is the uncertainty in- 
volved in the large x behavior of the deuterium data 22 : if this is determined as- 
suming a valence quark model 9 at large x, rather than fitting the data directly, the 
value 22 of r™ increases by a factor two; 9 however, the data are sufficiently uncer- 
tain that the two results are still compatible within errors. The problems related 
to small x are shown in Fig. 4, which displays the extrapolation of the E142 data, 
performed assuming a Regge-like behavior A™ ~ x 12 , 59 compared with an ex- 

x— »0 

trapolation of neutron data obtained by combining the SMC data 22 with the EMC 
proton data, 2 and then fitting a power behavior Eq. (3.8) to the last data point. 
These two extrapolations lead to values of T™ which differ by over 30%. 

The issue of scale dependence is particularly serious here because all neutron 
data are taken at very low Q 2 . This not only means that the (calculable) effects 
of QCD evolution due to difference in scale between x bins are now rather sizable 
(the correction on the first moment is of order 10%), 20 but also, that higher twist 
effects can be important. If these effects are included, then, for instance, the scale 
dependence of the nonsinglet first moment in Eq. (5.1) is modified as 



1 1 6 



CNs(a u - ad) + -^p- 



(5.4) 



The value of cht in Eq. (5.4) has been the subject of considerable controversy: 
QCD sum rules lead 63 to cht = —0.09 ±0.06, but a different estimate based on the 
same method 64 has cht = —0.15 ± 0.02; a bag model computation 65 even leads to 
a result with the opposite sign, cht = 0.16. The sum rule method appears to be 
self-consistent in that a different sum rule 66 leads essentially to the same result as 
Ref. 63; the theoretical uncertainty is perhaps of order 50%. 66 The scale dependence 
of the isotriplet Ti is displayed in Fig. 5; higher twist effects are included according 
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Ref. 


m 


Ti 


AE 


a s 


22 


4.7 


Tf = 0.023 ±0.025 


0.06 ±0.25 


-0.21 ±0.08 


59 


2.0 


ry = -0.022 ± 0.011 


0.57 ±0.11 


-0.01 ±0.06 


3 


5.0 


F[ = -0.055 ± 0.025 


0.24 ±0.23 


-0.11 ±0.08 


5 


3.0 


Tf = 0.044 ± 0.005 


0.35 ± 0.05 


-0.08 ±0.03 



Table 2: Summary of neutron and deuteron experimental results. All results hold at the given scale; 

r( = (rf + r?)/2 

to Eq. (5.4) with cut = —0.1. This shows that the magnitude of higher twist 
correction is comparable to that of three loop corrections, but smaller than the 
uncertainty on A. 

The neutron experiments are summarized in Table 2. The result of Ref. 3 
(shown in Fig. 4) is obtained by putting together the neutron data of Ref. 59 
with the values of g™ obtained combining g\ ± g™ determined from the deuterium 
experiment 59 with proton data. 2 ' 21 The values of AE and a s are as quoted by the 
respective references, except the entry corresponding to the preliminary result of 
Ref. 5, which we determined using the weak decay constants Eqs. (3.14)-(3.15), and 
the coefficient functions (3.17)-(3.18). 

Again, all experiments agree within errors. The error on the SMC data 22 
is mostly statistical, whereas the E142-E143 data, 59 ' 5 which have generally better 
statistics (but smaller kinematic coverage), have approximately equal statistical and 
systematic uncertainty. The significantly larger values of found by the SLAC 
experiments 59 ' 5 appear to be mostly due to the small-x region, as shown by the fact 
that the values of of Refs. 59 and 3 differ essentially because of the contribution 
from this region (notice however that they are compatible within errors). If the 
trend displayed by the SMC data at small x (Fig. 4) is confirmed, the small x 
extrapolation of the more precise data 5 would have to be corrected accordingly. 



6. Summary: the data and sum rules 

The determinations of Ti with proton, neutron, and deuterium targets can now 
be compared with each other and with QCD expectations; the results, evolved to 
a common scale using the perturbative coefficient functions of Sect. 3.2 to highest 
available order (omitting higher order estimates and higher twist corrections), are 
displayed in Fig. 6. 

First, the isotriplet first moment can be checked against the prediction of 
Eq. (5.1) (Bjorken sum rule 67 ), thereby checking isospin in this channel. As is 
apparent from Fig. 6, the prediction is perfectly verified within current experimen- 
tal errors, which means that it is verified to 10% accuracy. If, as the experimental 
precision improves, a violation of the sum rule within these bounds were found, it 
would not necessarily signal a breakdown of perturbative QCD: on the one hand, 
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it has been argued 68 that isospin breaking in meson-nucleon couplings could lead 
to a violation of the Bjorken sum rule up to 20%; also, the anomalous contribu- 
tion O [Eq. (4.3)] to the singlet matrix element AE, due to its infrared sensitivity, 
could develop 69 a nonsinglet component (triggered by SU(2) violation in light quark 
masses) of up to 10%. 

Then, a common value of AE can be extracted. This is best done 11 by fitting 
the expression (2.9) to the data with fixed values of the weak decay constants, and 
AE(Q 2 = oo) left as a free parameter. The result, obtained from all the data of 
Tables 1 and 2 (excluding the entries corresponding to Ref. 3) is 11 

AE(Q 2 = oo) = 0.31 ± 0.04. (6.1) 

This corresponds to a s (<5 2 = oo) = — 0.097 ± 0.018: the parton-based expectation 
that a s should be compatible with zero (Ellis- Jaffe sum rule 70 ) is thus violated by 
several standard deviations, with the theoretical implications discussed in Sect. 4.1. 

Finally, the considerable scale sensitivity of the first moments (compare Fig. 5), 
especially in the Q 2 < 5 GeV 2 region, can be used to measure a s . 71,11 The simplest 
way of doing this is to take advantage of the fact that both the scale dependence 
and the normalization of the isotriplet combination in Eq. (5.1) are accurately 
predicted, so that comparing an experimental determination of — r™ at finite 
Q 2 to the asymptotic value Eq. (3.14) gives immediately a determination of a s . 
Including all known perturbative corrections and estimates of higher loops (but not 
higher twist corrections), and using of Ref. 59 (as reanalysed in Ref. 8) and 
of Ref. 5 gives 71 

a s (M 2 ) = O.mjftgg, (6.2) 

in good agreement with the current world average, and with surprisingly small 
error. This error, however, does not include the uncertainty related to the higher 
twist correction: if it is taken into account using Eq. (5.4) with cht = —0.12 ± 
0.06 (according to the sum rule estimate discussed in Sect. 5), then 71 a s (M§) = 
O.IIS^o'q] 17 . A determination which includes all the available information can be 
arrived at fitting Eq. (2.9) to all the available data, with both AE and a s left as 
free parameters; the result is 11 (neglecting again higher twist effects) 

AE(Q 2 = oo) = 0.33 ±0.04; a s {M 2 z ) = 0.125 ± 0.006. (6.3) 

This result is once more surprisingly accurate, and in good agreement with Eq. (6.2). 
It would be interesting to also fit simultaneously a s , AE, and the octet combination 
+ — 2a s , but current data are not accurate enough. 71 ' 11 
In conclusion, combining all available results leads to a consistent picture, in 
excellent agreement with the predictions of perturbative QCD, and with interesting 
implications for the nucleon structure. While the improving quality of experimental 
data requires a more sophisticated theoretical analysis, it is now possible to test our 
understanding of several subtle perturbative QCD effects, and to start coping with 
the nonperturbative effects which determine the structure of the nucleon. 
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Figure Captions 



Fig. 1: Experimental determinations of AE P (at the scale of the respective experi- 
ments, see Tab.l below): a) cross, EMC (1988) ^circle, EMC (1989); 2 diamond, 
world average (1994); 3 square, SMC (1994); 4 star, E142 (1994). 5 ' 6 b) Theoret- 
ical reanalysis of the 1989 experiment 2 : cross, published value 2 ; circle, Jaffe 
and Manohar (1990); 7 diamond, Ellis and Karliner (1993); 8 square, Close and 
Roberts (1993). 9 

Fig. 2: (from Ref. 20): a) Scale dependence of the asymmetry Ai(x, Q 2 ), compared to 
the data. The crosses indicate data from Ref. 21 and the squares data from 
Ref. 2; the low (high) curves correspond to the lower (upper) edge of the x-bin, 
the solid (dashed) lines corresponds to a maximal polarized gluon (no polarized 
gluon). b) Effect on the data 21 ' 2 of the scale dependence correction. 

Fig. 3: a) (from Ref. 2) Experimental data for g{(x) at (Q 2 ) = 10.7 GeV 2 (filled dots, 
EMC data; 2 open circles, SLAC data 21 ); the dashed line is the small-x extrap- 
olation. b)(from Ref. 4) The data of Fig3a (here as triangles) compared to the 
new data. 4 

Fig. 4: (from Ref. 3) The E142 59 and SMC/EMC 4 ' 2 data with the respective small x 
extrapolations (dashed: E142; solid: SMC/EMC) 

Fig. 5: (adapted from Ref. 20) Scale dependence of the isotriplet first moment 
Eq. (5.1). Dashdot: tree-level result (Cns = 1); dotted: Leading order; dashed: 
three loops; solid: with higher twist. The three sets of curves correspond to 
the values A = 267, 383, 509 MeV. 

Fig. 6: Comparison of experimental determinations of F± for proton, 5 neutron 59 and 
(dashed) deuteron. 5 The dotted lines indicate alternative proton 2 and neutron 3 
determinations. The Bjorken sum rule prediction Eq. (5.1) is also shown (dot- 
dashed). All results have been evolved to Q 2 = 5 GeV 2 . 
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